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Abstract 

Anomalies are known to appear in the perturbation theory for the one-dimensional 
Anderson model. A systematic approach to anomalies at critical points of products of 
random matrices is developed, classifying and analysing their possible types. The associ- 
ated invariant measure is calculated formally. For an anomaly of so-called second degree, 
it is given by the groundstate of a certain Fokker-Planck equation on the unit circle. The 
Lyapunov exponent is calculated to lowest order in perturbation theory with rigorous 
\ control of the error terms. 
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1 Introduction 



Anomalies in the perturbative calculation of the Lyapunov exponent and the density of states 
were first found and analysed by Kappus and Wegner |KWj when they studied the center of 
the band in a one- dimensional Anderson model. Further anomalies, albeit in higher order 
perturbation theory, were then treated by Derrida and Gardner |DGj as well as Bovier and 
^ ■ Klein [BK . More recently, anomalies also appeared in the study of random polymer models 
jJSSj . Quite some effort has been made to understand anomalies in the particular case of the 
Anderson model also from a more mathmatical point of view |CK|ISVW] . However, Campanino 
and Klein |CKj need to suppose decay estimates on the charcacteristic function of the random 
potential, and Shubin, Vakilian and Wolff |SVWj appeal to rather complicated techniques from 
harmonic analysis (allowing only to give the correct scaling of the Lyapunov exponent, but not 
a precise perturbative formula for it). 

It is the purpose of this work to present a more conceptual approach to anomalies of products 
of random matrices. In fact, various types may appear and only those of second degree (in the 
sense of the definition below) seem to have been studied previously. Indeed, this is the most 
difficult and interesting case to analyse, and the main insight of the present work is to exhibit 
an associated Fokker-Planck operator, the spectral gap of which is ultimately responsible for 
the positivity of the Lyapunov exponent. In the special case of the Anderson model, a related 
operator already appeared in BKj . Here it is, however, possible to circumvent the spectral 
analysis of the Fokker-Planck operator and prove the asymptotics of the Lyapunov exponent 
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more directly ( cf. Section 15. 4|) . The other cases of various first degree anomalies are more 
elementary to analyse. Examples for different types of anomalies are given in Section |3 

Acknowledgement: This is a contribution to the conference procedeings of OTAMP2004 held 
in Bedlewo in June 2004. First of all, I would like to thank the organizers for the invitation 
to the wonderful and inspiring conference in the convivial atmosphere of Bedlewo, allowing for 
quite unexpected encounters; second of all, many thanks to my new collegues Andreas Knauf 
and Hajo Leschke for a few patient discussions on the matters of the present work. 

2 Definition of anomalies 

Let us consider families (T A CT ) AeIR)0 - e2 of matrices in SL(2,R) = {T e Mat 2X 2(R) I det(T) = 1} 
depending on a random variable a in some probability space (£, p) as well as a real coupling 
parameter A. In order to avoid technicalities, we suppose that p has compact support. The 
dependence on A is supposed to be smooth. The expectation value w.r.t. p will be denoted by 
E. 

Definition 1 The value A = is anomaly of first order of the family (T\ iCr )AeK,o-es if for all 
a E S; 

T , CT = ±1 , (1) 

with a sign that may depend on o e S. In order to further classify the anomalies and for later 
use, let us introduce P a ,Qa £ sl(2,R) by 

MT\^M~ l = ±exp(AP CT + A 2 Q a + 0(X 3 )) , (2) 

where M e SL(2, R) is a A- and a -independent basis change to be chosen later. An anomaly is 
said to be of first degree ifE(P a ) is non-vanishing, and then it is called elliptic z/det(E(P cr )) > 0, 
hyperbolic if det(E(P CT )) < and parabolic if det(E(P CT )) = 0. Note that all these notions are 
independent of the choice of M. 

If E(P cr ) = 0, but the variance of P a is non-vanishing, then an anomaly is said to be of 
second degree. 

Furthermore, for k e N, set a = (o~(k), . . . , c(l)) G S = S xfc , as well as p = p xfe and 
Pa.ct = Tx,er{k) " " " Tx,er(i). Then A = is anomaly of kth order of the family (Tx t(T )\eM.,aes if the 
family (T\^) XeU g. e g has an anomaly of first order at A = in the above sense. The definitions 
of degree and nature transpose to kth order anomalies. 

As is suggested in the definition and will be further explained below, we may (and will) 
restrict ourselves to the analysis of anomalies of first order. In the examples, however, anomalies 
of higher order do appear and can then be studied by the present techniques (cf. Section^. 
Furthermore, by a change of variables in A, anomalies of degree higher than 2 can be analysed 
like an anomaly of second degree. 
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Anomalies are particular cases of so-called critical points studied in |JSS[ ISSSj . namely 
A = is by definition a critical point of the family (T\ t(T ) AeK.o-eE if for all a, o' G E: 

[To^Toy] = 0, and |Tr(T 0j(T )| < 2 or T , CT = ± 1 . (3) 

Critical points appear in many applications like the Anderson model and the random polymer 
model. In these situations anomalies appear for special values of the parameters, such as the 
energy or the coupling constant, cf. Section 



3 Phase shift dynamics 

The bijective action St of a matrix T G SL(2, R) on S 1 = [0, 2tt) is given by 

This defines a group action, namely Stt' = StSt>- In particular the map St is invertible and 
S^ 1 = St-i- Note that this is actually an action on KP(1), and S 1 appears as a double cover 
here. In order to shorten notations, we write 



S\,<t — S MTx a M 



-i 



Next we need to iterate this dynamics. Associated to a given semi-infinite code uo = (o" n )n>i 
with cr n G E is a sequence of matrices (Tx jan )n>i- Codes are random and chosen independently 
according to the product law p (XlN . Averaging w.r.t. p® n is also denoted by E. Then one defines 
iteratively for iV G N 

= SxrAS&W . S U°) = e - (5) 

This is a discrete time random dynamical system on S 1 . Let us note that at an anomaly of first 
order, one has S\ ta (6) = 9 + 0(X) or S\ >a (9) = 9 + n + 0(\) depending on the sign in (JTJ). As all 
the functions appearing below will be 7r-periodic we can neglect the summand 7r, meaning that 
we may suppose that there is a sign + in JTJ for all a (this reflects that the action is actually 
on projective space). 

In order to do perturbation theory in A, we need some notations. Introducing the unit 

V2 [ - % 

p a {6) = 3m ( N " r ) , Q<t{9) = 3m 



vector v = 4= ( ^ ) , we define the first order polynomials in e 2lB 
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as well as 

a a = (v\P a \v) , f} a = (v\P a \v) 
Hence p a (9) = ^sm(a a — (3 a e 2%e ). Now starting from the identity 
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(v\MT a M- 1 \e ) ' 

the definition (J2J) and the identity (v\e$) = -j= e ld , one can verify that 



V. (v|ee) 2 (v\e ) 2 (w|e e ) 2 / 

As one readily verifies that 

P 2 = -det(P CT ) 1 , 3m ( {V ] Kle f ) = - P(r (9)dep a (9) , 

it follows that 

Sx,M = + A P(7 (0) + \ 2 q„(6) + i \ 2 p a (9) d e pa(0) + 0(A 3 ) . (6) 
Finally let us note that 

550) = * - A P<r ((9) - A 2 ^) + iAV(0)&M0) + 0(X 3 ) , (7) 

as one verifies immediately because S\ t(T (S x 1(9)) = $ + ^(A 3 ), or can deduce directly just as 
above from the identity exp (A P a + A 2 Q a + ^(A 3 ))^ 1 = exp (-A P a -\ 2 Q a + C(A 3 )). 



4 Formal perturbative formula for the invariant measure 

For each A, the family (MT x ^M~ l ) ae Y, and the probability p define an invariant probability 
measure v x on S 1 by the equation 

J du x (9) f(9) = E J du x (9) f(SxA0)) , / e C(S l ) . (8) 

Furstenberg proved that this invariant measure is unique whenever the Lyapunov exponent of 
the associated product of random matrices (discussed below) is positive (e.g. |BLj ) and in this 
situation v\ is also known to be Holder continuous, so, in particular, it does not contain a point 
component. For the study of the invariant measure at an anomaly of order k, it is convenient 
to iterate (jSJ): 

J du x (9) f(6) = E J du x (9) f(S k x J9)) , / e CiS 1 ) . 

Replacing (E, p) by (£, p) therefore shows that the families (T AjCr ) Ag]K)(TgS and (T XjS .) xm . j, have 
the same invariant measure. Hence it is sufficient to study anomalies of first order. 



4 



The aim of this section is to present a formal perturbative expansion of the invariant measure 
under the hypothesis that it is absolutely continuous, that is dv x (9) = pxiff)^- with p x = 
po + Xpi + 0{\ 2 ). Then © leads to 

E (d e S^{9)p x {S^{9))) = Px (9), (9) 

which with equation (J7J) gives 

p x - Xd e (E(p a )p^ + \ 2 ^d e (E(pl)d ePx + E( Pa d ePl7 ) p x - 2E( g(T )p A ) + G(A 3 ) = Px . 

We first consider an anomaly of first degree. As E(P CT ) 7^ 0, it follows that E(p CT ) is not 
vanishing identically. Therefore the above perturbative equation is non-trivial to first order in A, 
hence F>( Pa )p should be constant. If now the first degree anomaly is elliptic, then det E(P a ) > 
which can easily be seen to be equivalent to |E(a CT )| > |E(/3 CT )|, which in turn is equivalent to 
the fact that E(p a (9)) does not vanish for any 9 G S 1 . For an elliptic anomaly of first degree, 
the lowest order of the invariant measure is therefore 

P ° = E^)' 

with an adequate normalization constant ceR. If on the other hand, the anomaly is hyperbolic 
(resp. parabolic), then E( Prr ) has four (resp. two) zeros on S 1 . In this situation, the only 
possible (formal) solution is that po is given by Dirac peaks on these zeros (which is, of course, 
only formal because the invariant measure is known to be Holder continuous). In Section I5~21 
we shall see that for the calculation of certain expectation values w.r.t. the invariant measure, 
it looks as if it were given by a sum of two Dirac peaks, concentrated on the stable fixed points 
of the averaged phase shift dynamics. These fixed points are two of the zeros of E( PrT ). 

Next we consider an anomaly of second degree. As then E(po-) = 0, it follows that the 
equation for the lowest order of the invariant measure is 

- d e (e( P 1) d eP o + EfadopJ po - 2 E(q a ) p ) = . 

Now E(p^.) > unless p-almost all p a vanish simultaneously for some 9, a (rare) situation which 
is excluded throughout the present work. Then this is an analytic Fokker-Planck equation on 
the unit circle and it can be written as L po = where £ is by definition the Fokker-Planck 
operator. Its spectrum contains the simple eigenvalue with eigenvector given by the (lowest 
order of the) invariant measure p calculated next. Indeed, 

i E(pl) d e po + ^ E(p a de P a) Po ~ E{q a ) p = C , 

where the real constant C has to be chosen such that the equation admits a positive, 2n- 
periodic and normalized solution p - It is a routine calculation to determine the solution using 
the method of variation of the constants. Setting 
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it is given by 



-k(2tt) 



if(27r) 



r p «(0) 

^ = ^2^ l (C, ^ ) + 1) ' (10) 
E(p2(0))a 

where c is a normalization constant. It is important to note at this point that po($) is an 
analytic function of 9. 

The rest of the spectrum of C is discrete (C has a compact resolvent), at most twice degen- 
erate and has a strictly negative real part, all facts that can be proven as indicated in [RisJ. As 
already stated in the introduction, we do not need to use this spectral information directly. 



5 The Lyapunov exponent 

The asymptotic behavior of the products of the random sequence of matrices (T\^ n ) n >i is 
characterized by the Lyapunov exponent |BL[ A. III. 3. 4] 



7 (A) = lim — E log 

N— >oo iv 



N 



n=l 



(11) 



where 9 is an arbitrary initial condition. One may also average over 9 w.r.t. an arbitrary 
continuous measure before taking the limit [JSS, Lemma 3]. A result of Furstenberg states a 
criterion for having a positive Lyapunov exponent [BLJ. A quantitative control of the Lyapunov 
exponent in the vicinity of a critical point is given in Proposition 1], however, only in the 
case where the critical point is not an anomaly of first or second order. The latter two cases 
are dealt with in the present work. 

Let us first suppose that the anomaly is of first order. Because the boundary terms vanish 
in the limit, it is possible to use the matrices MT A o . n M _1 instead of T Aj(Tn in ffTTj) . Furthermore, 
the random dynamical system (jSJ) allows to expand (JTTJl into a telescopic sum: 



7(A) 



j N-l 

J™ N £ E l0 S (|| MT W 1 M- 1 e 5 ^ (e) ||) 

n=0 



(12) 



Up to terms of order 0(X 3 ), we can expand each contribution of this Birkhoff sum: 



\og(\\MT Xta M~ l e e \\) = A (e e \P a \e e ) + l - A 2 ((e e | (|P CT | 2 + Q a + P 2 a )\e e ) - 2 {e e \P a \e e ) 

= l -^e [2\{3 a e™ + A 2 (|&| 2 + (v\ (\P a \ 2 + 2 Q a ) \v)e™ - /3 2 e 4 ^) 
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where we used the identity 

(e e \T\eg) = ^ Tr(T) + {(v\T\v) e™) , 

holding for any real matrix T, as well as Tr(P CT ) = Tr(Q a ) = and Tr(|P CT | 2 + P 2 ) = 4 \f3 a \ 2 . 
Let us set 

i 

UN) = ± £ E(e 2 ^)) , j = 1,2, (13) 

n=0 

and introduce Jj by Ij(N) = Ij + 0(A) for iV sufficiently large. We therefore obtain for an 
anomaly of first order 

7 (A) = ^E^e(2\p a h + \ 2 {\p a \ 2 + (v\(\P <T \ 2 + 2Q <T )\v)h - #I 2 )) + 0(A 3 ) . (14) 

For an anomaly of second order, one regroups the contributions pairwise as in Definition Q 
namely works with the family (T x ,a)xm,9e^ 2 where T\ t $ = Ta,<t(2)?V(i) for a = (a (2), a(l)), 
furnished with the probability measure p = p x2 . This family has an anomaly of first order, 
and its Lyapunov exponent is exactly twice that of the initial family (Tx^xm^es- It is hence 
sufficient to study anomalies of first order. 

5.1 Elliptic first degree anomaly 

Let us consider the matrix E(d\T\ t(T \\ =Q ) Gsl(2,M). For an elliptic anomaly, the determinant 
of this matrix is positive. Its eigenvalues are therefore a complex conjugate pair ±z |, so that 
there exists a basis change M G SL(2,R) such that 

E(M9 A T A , CT | A=0 M- 1 ) = E(P CT ) = i ( ° ) . (15) 
It follows that E(a CT ) = z f and E(/3 CT ) = 0, implying 7(A) = 0(X 2 ). Furthermore, from (jSJ, 
e awSA,aW = (1 + 2zj A 3m(a, - Z^e 3 *)) e 2 ^ + C(A 2 ) . 

Hence 

1 

n=0 

so that /j(iV) = 0(X, (NX)- 1 ) and I, = 0. Replacing into (JUJ), this leads to: 
Proposition 1 If X = is an elliptic anomaly of first order and first degree, then 

7(A) = I A 2 E(|/3 CT | 2 ) + 0(A 3 ). 
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In order to calculate f3 a in an application, one first has to determine the basis change (J15j) . 
then P a before deducing (3 a and Ed/^l 2 ), see Section|H]for two examples. Let us note that after 
the basis change, (JT5j) implies E(p CT (0)) = r]. Hence the invariant measure is to lowest order 
given by the Lebesgue measure after the basis change. 

The term 'elliptic' indicates that the mean dynamics at the anomaly is to lowest order a 
rotation. For a hyperbolic anomaly it is an expansion in a given direction and a contraction 
into another one. These directions {i.e. angles) can be chosen to our convenience through the 
basis change M, as will be done next. 

5.2 Hyperbolic first degree anomaly 

For a hyperbolic anomaly of first degree, the eigenvalues of E(<9aTa j(7 |a=o) are ±^ and there 
exists M e SL(2, R) such that 

E(M9 A T A , (T | A=0 M- 1 ) = E(P CT ) = 5 ( {J ^ ) • (16) 
It follows that E(a a ) = and E(/3 ff ) = |, so that dHJ) leads to 

7 (A) = \ \iiMh) + 0(X 2 ). 

We now need to evaluate I\. Introducing the reference dynamics S x {9) = 9 — |A/isin(2#) as 
well as the centered perturbation r a {6) = ?sm(a a — (/3 CT — |) e 2ld ), it follows from © that the 
phase shift dynamics is 

Sx,*(8) = Sx(6) + Xr a (9) + 0(X 2 ) . (17) 

The non-random dynamics S\ has four fixed points, 9 = 0, |, tt, If A/x > 0, 6 1 = and 7r are 
stable, while = | and ^ are unstable. For A/i < 0, the roles are exchanged, and this case 
will not be considered here. Unless the initial condition is the unstable fixed point, one has 
S^(9) — > 0, 7r as n — > oo. Furthermore, if 9 is not within a (9(A2)-neighborhood of an unstable 
fixed point, it takes n = C(A~2) iterations of S\ in order to attain a C(A2)-neighborhood of 0. 
We also need to expand iterations of S\: 

S k x (9 + Xr a (9') + 0(X 2 )) = S k x (9) + Xd e S k x (9) r a (9') + 0{X 2 ) , 

where the corrective term 0(X 2 ) on the r.h.s. is bounded uniformly in k as one readily realizes 
when thinking of the dynamics induced by S\. Furthermore, dgS x (9) = 0(1) uniformly in k. 
Iteration thus shows: 

n 

SIM = S n x (9) + A^^r^^^- 1 ^)))^^- 1 ^)) + 0(nX 2 ). 

k=l 
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Let us denote the coefficient in the sum over k by s&. Then Sk is a random variable that depends 
only on a% for / < k. Moreover, Sk is centered when a conditional expectation over Ok is taken. 
Taking successively conditional expectations thus shows 

E(e 2l V)) = e 2 ^ (E(e 24A ^^ lsft + °( nA2 )) + 0(A 2 )) = e 2 ^"W + 0(nA 2 ) . (18) 

3 

Choosing n = X~2 gives according to the above 

E(e 2l5 U e) ) = 1 + , (19) 

unless # is within a (9(A5)-neighborhood of | or In the latter cases, an elementary argument 
based on the central limit theorem shows that it takes of order E(r (T ( : |) 2 )A~i iterations to diffuse 
out of these regions left out before. Supposing that one does not have r a (~) = for p-almost 
all a, one can conclude that (fT9"|) holds for all initial conditions 9. Consequently I\ = 1 + 0(X^) 
so that: 

Proposition 2 If X = is an hyperbolic anomaly of first order and first degree, andr a (j) does 
not vanish for p-almost all a, one has 

7 (A) = 1|Aai| + O(Xl). 

The argument above shows that the random phase dynamics is such that the angles S^^iO) 

are for most n and u in a neighborhood of size A^ of the stable fixed points 9 = 0, it. This does 
not mean that for some n and to, the angles are elsewhere; in particular, the rotation number 
of the dynamics does not vanish. However, this leads to corrections which do not enter into 
the lowest order term for the Lyapunov exponent. 

5.3 Parabolic first degree anomaly 

This may seem like a pathological and exceptional case. It turns out to be the mathematically 
most interesting anomaly of first degree, though, and its analysis is similar to that of the 
Lyapunov exponent near band edges (this will be discussed elsewhere). First of all, at a 
parabolic anomaly of first degree, there exists M eSL(2,R) allowing to attain the Jordan 
normal form: 

E(Md x T x<a \ x=0 M- 1 ) = E(P CT ) = f o ) * (20) 

Thus E(a CT ) = E(/3 cr ) = — 1 1 so that 7(A) = ^XQm(Ii) + 0(X 2 ). Introducing the reference 
dynamics <S\{0) = 9 + |(cos(2#) — 1) as we U as the centered perturbation r a {9) = Qm(a a + 
\ — (Pa + f ) e 2i61 ), the dynamics can then be decomposed as in (|T7|) . Moreover, the argument 
leading to (|18j) directly transposes to the present case. However, this does not allow to calculate 

3 

the leading order contribution, but shows that 7(A) = O^X?). 
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5.4 Second degree anomaly 



At a second degree anomaly one has E(/3 CT ) = 0, so in order to calculate the lowest order of the 
Lyapunov exponent one needs to, according to (JHJ), evaluate l\ and J 2 . For this purpose let us 
introduce an analytic change of variables Z : S 1 — > S 1 using the density p given in (jTUj) : 

9 = Z{9) = [ dB' p {6') . 
Jo 

According to Section |U one expects that the distribution of 9 is the Lebesgue measure. We 
will only need to prove that this holds perturbatively in a weak sense when integrating analytic 
functions. 

We need to study the transformed dynamics S\ jCT = Z o S\. a o Z^ 1 and write it again in the 
form: 

Sx, a {§) = 9 + Xp a (9) + \ 2 q a {9) + ±\ 2 p a (9)d § p a (9) + G(A 3 ) . (21) 
As, up to order C(A 3 ), 



ZoS^(9) = Z(9) + (\vM + \\ 2 V*{0)d e pM + X 2 q a (9^d e Z(9) + ^ X 2 p 2 a d 2 e Z(9) 



one deduces from 

2 / 1 

d e Z = p , d 2 e Z = — — (C - - E(p a d e p a ) p + E(q a ) p 
that, with 9 = Z- l {9), 

p a (9) = Pa (9)p (9) , q„(9) = q„(9)p (9) . 
A short calculation shows that the expectation values satisfy 

E(p a (9)) = , E(q a (9)) = I E (p a {9) d § p a (9)) - C , (22) 

where C is as in (fTUj) . 

Given any analytic function f on S 1 , let us introduce its Fourier coefficients 

trnO 



p2tt m 



There exist a, £ > such that f m <ae ^' m L We are interested in 



AT-l 



j/W 4 E E Mo. 



AT 

n=0 
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where, for sake of notational simplicity, we introduced 6 n = S^^ip) for iterations defined just 
as in (J^J). 

Lemma 1 Suppose E(p 2 ) > and that f is analytic. Then 

I f <N) = f + 0(X, (X'N)- 1 ) , 
with an error that depends on f . 

Proof. Set f = ~E(p 2 ). This is an analytic function which is strictly positive on S 1 . Further- 
more let F be an auxilliary analytic function with Fourier coefficients F m . Then, using (|21j) 
and the identities (|22p. 

1 N ~ 1 

I p (N) = — E ^ ^ F m e 4m4 (l - zmCA 2 + mX 2 {id § - m)r(6 n ) + C(m 3 A 3 , iV" 1 )) 

n=0 meZ 

= i p (N) + \ 2 L c p, + ^, +f p,(N) + OiX^N- 1 ) , 
where the prime denotes the derivative. Hence we deduce that for any analytic function F 

L c p +if p, y (N) = 0(X, (A 2 A^) _1 ) . (23) 

Now, extracting the constant term, it is clearly sufficient to show I AN) = 0(X, (A 2 A^) _1 ) 
for an analytic function / with f = 0. But for such an / one can solve the equation 

/ = - CF' + (fFj (24) 



for an analytic and periodic function F and then conclude due to ([23)1 . Indeed, by the method 
of variation of constants one can always solve (|24|) for an analytic F' . This then has an 
antiderivative F as long as F' does not have a constant term, i.e. the zeroth order Fourier 
coefficient of the solution of (J24j) vanishes. Integrating (J24j) w.r.t. 6, one sees that this is 
precisely the case when f Q = as long as C ^ 0. If on the other hand (7 = 0, then (J24)) can 
be integrated once, and the antiderivative J f of / chosen such that f _1 J f does not have a 
constant term. Then a second antiderivative can be taken, giving the desired function F in this 
case. □ 

In order to use this result for the evaluation of Ij(N) defined in (I13|) . let us note that 



j 

i ^-i 

n=0 



e 2ijz-\e n ) 



Hence up to corrections the result is given by the zeroth order Fourier coefficient of the analytic 
function f(6) = e 2t ^ z so that after a change of variables one gets: 

/■27T in 

I S (N) = / - Po (0)e^ e + 0(X, (X 2 Ny 1 ) . 
Jo ^ 



11 



Proposition 3 If A = is an anomaly of first order and second degree and E(p 2 ) > 0, then 
one has, with p given by ([TU]). 

7(A) = Y^C^^ + Vm(\P«\ 2 + 2QM)e™-E(tt)e^] + 0(A 3 ). 

In the above, anomalies of first degree were classified into elliptic, hyperbolic and parabolic. 
Second degree anomalies should be called 'diffusive' (strictly if E(p 2 ) > 0). The random 
dynamics of the phases is diffusive on S 1 , with a varying diffusion coefficient and furthermore 
submitted to a mean drift, also varying with the position. It does not seem possible to transform 
this complex situation into a simple normal form by an adequate basis change M. 

Let us note that the Lyapunov exponent at an anomaly does depend on the higher order 
term Q a in the expansion (J2J), while away from an anomaly it does not depend on Q a [SEE] . Of 
course, the coefficient of A 2 in Proposition|3]cannot be negative. Up to now, no general argument 
could be found showing this directly (a problem that was solved in jSSSl Proposition 1] away 
from anomalies). For this purpose, it might be of help to choose an adequate basis change M. 

6 Examples 

6.1 Center of band of the Anderson model 

The transfer matrices of the Anderson model are given by 

rv=(V" D '), (25) 

where v a is a real random variable and E e M is the energy. The band center is given at E — 0. 
In order to study the behavior of the Lyapunov exponent at its vicinity, we set E = eA 2 for 
some fixed e G M. Then the associated family of i.i.d. random matrices has an anomaly of 
second order because 

T„ = = - exp(A( «*■ ) + A2 (-!>WV, i ) + 0(V) ) , 

where a — (a 2 ,cr 1 ). It follows that a& = J-(w<j 2 + v ai ) and (3„ = 2~( v °-2 ~ v ^i)- ^ now v <? is n °t 
centered, then one has an elliptic anomaly of first degree and Proposition Q combined with the 
factor 1/2 due to the order of the anomaly, implies directly (no basis change needed here) that 

7(A) = ^ A 2 (E( V 2 ) - EK) 2 ) + 0(X 3 ). 

If, on the other hand, v a is centered, one has a second degree anomaly and can apply Propo- 
sition El One readily verifies that E(|/?5-| 2 ) = |E(t> 2 ) and E(/3f) = — |E(t> 2 ), and furthermore 
that the second term in Proposition El always vanishes, so that 
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i i-2-n in 

7 (A) = - A 2 E(vl) / - p (6) (1 + cos(4£)) + O(X') , 

which is strictly positive unless v a vanishes identically (it was already supposed to be centered). 
If one wants to go further, E(|pa-(#)| 2 ) = |E(t> 2 )(l + cos 2 (2#)). In the case e = 0, one then 
has E(Q^) — so that po(9) = c(l + cos 2 (2#))~^ with a normalization constant c that can be 
calculated by a contour integration. This proves the formula given in |DGj . 



6.2 A particular random dimer model 

In the random dimer model, the transfer matrix is given by the square of (|25jl . with a potential 
that can only take two values Xv a = crv where i)6l and o G { — 1, 1} (e.g. [JSSJ). Hence for a 
given energy E G K it is 



av-E -1 x 2 



T = 

Now the energy is chosen to be E = v + A. Then 7\ )Cr = T^ +A has a critical point if |u| < 2. 
For the particular value v = fixing hence a special type of dimer model, one has: 

This family has now an anomaly of second order and first degree because 

One readily verifies that the determinant of E(M _1 P (T M) is equal to 7 — 2E(cr) — E(cr) 2 and 
hence positive so that the anomaly is elliptic. Therefore the general result of Section 15.11 can 
be applied. Let us set e = E(er). The adequate basis change (without normalization of the 
determinant) is 



M 



V7-2e- e 2 
\/2(l-e) 3-e 



A calculation then gives 



1 / 2^2(^-6) (a-3) v / 7-2e-e 2 



4(l- e )( -e)-(e-3)(7- -e-e 2 ) -2x/2((T 



allowing to extract (3 a and then E(|/3o-| 2 ), leading to (this contains a factor 1/2 because the 
anomaly is of second order) 



2 2(l-e 2 ) (. 2(e-l« 
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Note that if e = 1 or e = — 1 so that there is no randomness, the coefficient vanishes. This 
special case was left out in [JS5] (the condition E(e 4ir? ) 7^ 1 in the theorems of jJSSj is violated). 
Within the wide class of polymer models discussed in |JSSj . models with all types of anomalies 
can be constructed, and then be analyzed by the techniques of the present work. 
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